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ABSTRACT 

If a component of cosmological dark matter is made up of massive particles - 
such as sterile neutrinos - that decay with cosmological lifetime to emit photons, 
the reionization history of the universe would be affected, and cosmic microwave 
background anisotropies can be used to constrain such a decaying particle model 
of dark matter. The optical depth depends rather sensitively on the decaying 
dark matter particle mass mam, lifetime Td m , and the mass fraction of cold dark 
matter / that they account for in this model. Assuming that there are no other 
sources of reionization and using the WMAP 7-year data, we find that 250 eV 
< m dm < 1 MeV, whereas 2.23 x 10 3 yr < r dm // < 1.23 x 10 18 yr. The best fit 
values for m dm and r dm // are 17.3 keV and 2.03 x 10 16 yr respectively. 

Subject headings: dark matter, sterile neutrinos 



-3 - 



Introduction 



There have been tremendous progress in cosmology in the past decade. The avai lability 



20101 ) has 



of high quality observational data such as those from WMAP (IKomatsu et al. 
led to tight constraints on cosmological parameters and models. There is now a standard 
model of cosmology, in which only a small portion of the total mass-energy in the universe 
is ordinary matter, the rest being dark components which we have litt le understanding of. 



In so me of the proposed dark matter models, such as sterile neu trino flD odelson & Wi drow 



Cen 



1994) , the dark matter particles may decay and emit photons ( IBorzumati et al. 



2008 



20011 ). which are redshifted with the expansion of the universe and may eventually 



ionize hydrogen and helium at later times. Therefore, decaying dark matter particles may 
contribute to reionization and imprint their signatures on the cosmic microwave background 
anisotropies (CMBA). In this paper, we constrain the mass and life time of decaying dark 
matter particles by using the WMAP data of CMBA. 



There are strong evidences for reionization at late universe ( 



sources of which have been pro posed, such as star 



UV radiation from bl a ck holes ( Sasaki fc Umemura 



( iTegmark et al. 



1993|). 



Becker et al. 



200 



ormat ion (IGnedin fc Ostriker 



), m any 



19971 ) 



19961 ). and supernova- driven winds 



Biermann fc Kusenkol (120061 ) point out that the X-ray photons 



produced in the decays of sterile neutrinos can boost the production of molecular hydrogen, 
and as a result the rates of cooling of gas and early star form ation are increas e d, lea ding 



to reionization at redshift consistent with the WMAP results. 



Boyarsky et al. 



(120061 ) use 



extragalactic diffuse X-ray ba ckground to constra in the decay rate of sterile neutrinos as a 



warm dark matter candidate. 



Seljak et al. 



(120061 ) use the Ly-alpha forest power spectrum 



measured by the Sloan Digital Sky Survey and high-resolution spectroscopy observations 
in combination with WMAP data and galaxy clustering to constrain sterile neutrino 
masses. The lower limits obtained are 13.1 keV at 95% CL. and 9.0 keV at 99.9% C.L.. 
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In IZhang et al.l (120071 ) . decaying dark matter is also considered to be an energy source 
of reionization. However several approximations are made in that paper; in particular, 
the fraction of the decay energy deposited in baryonic gas is simply characterized by 
a phenomenological parameter. In this paper, we do not make approximations about 
the amount of energy absorbed by the baryons. The ionization and heating rates are 
calculated using the appropriate cross sections. Furthermore, we vary both the decaying 
dark ma tter particle param eters and cosmological parameters to fit the CMBA data, 



while in 



Zhang et al 



( 120071 ) only the decaying dark matter particle parameters and the 



scalar amplitude are vari ed to fit the CMBA spectrum. There is also another earlier work 



( jMapelli fc Ferrara 



20051 ) studying the effect on reionization by the decaying process of 
sterile neutrinos, where a relation between sterile neutrino mass and lifetime is used, which 
is based on the assumption that sterile neutrinos are the dominant component of dark 
matter. However in this paper we do not assume any relation between the decaying sterile 
neutrino mass and lifetime □; we treat them as two independent parameters instead, and we 
introduce another free parameter, the mass fraction of dark matter that are decaying, /. 

The optical depth depends rather sensitively on the mass md m and life time Td m of dark 
matter particles, as well as /. However, to good approximation, the effects of / and Td m 
are degenerate and only their ratio is an independent parameter. We then constrain these 
parameters, md m and r dm // by the WMAP 7- year data. Assuming that such a decaying 
process is the only source of reionization [^|, we find that m^ m is less than about 1 MeV, and 



x The lifetime refers to the radiative channel only. Sin ce the sterile neutr inos can also 



19951 ) . 



decay into 3 active neutrinos, the total lifetime ~ Td m /128 (IBarger et al. 

2 There is a well-known discrepancy between reionization redshifts deduced from CMBA 
and quasar absorption line observations. However, the constraints based on quasar absorp- 



tion line observations are highly model-dependent (IMcGreer et al 



201ll ). In particular, the 
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T dm/ f is less than about 10 18 yr, with the best-fit values being 17.3 keV and 2.03 x 10 16 yr 
respectively. 

In Section 2, we present the calculation of the ionization fraction and optical depth 
in this decaying dark matter model. The Markov Chain Monte Carlo fitting to WMAP 
data results and discussion are presented in Section 3, and Section 4 is a summary and 
conclusion. 



2. The Model 



The evolution of the ionization fractions xk(z) and x^ e (z) (for hydrogen and helium 
respectively) and the matter temperature T m (z) satisfy the coupled ordinary differential 
equations 

dxft(z) —1 



dz 



dx-uJz 



[1 + z)n u (z)H(z 
-1 



-[R m (z) + R sxH (z)] 

,. . , W7 , ARiUc(z) + RsxRe(z)] 

dz (1 + z)nu c (z)H(z) 

= -2 Rhujz) + RhKejz) + R s t(z) 

dz 3fc B (l + z)H(z) [1 + x H (z) + xn c (z)}n R (z) 



dT m (z) 



(1) 
(2) 
(3) 



steep rise in the Gunn-Peterson effective optical depth at z > 6 is hig 



ily controversial, a s 



it is very sensitive to the assumed density field and continuum fitting ( IBecker et al. 



20071 ). 



Recently, direct and model-independent limits on the fraction of neutral Hydrogen at z w 5-6 
were obtained using the simple statistic of the covering fraction of dark pixel s, and they can 



be co nsistent with the ionization history derived from CMBA observations (IMcGreer et al. 
201ll ). On the other hand, a recent work shows that model independent joint CMBA- 
quasar absorption line constraints still permit a broad range of reionization history for z > 6 



(IMitra et al. 



20121 ) . 
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where R sx {z) and R s t(z) are the standard net recombination and net heating rates 
respectively, and Ris(z) and Rhs(z) are the additional ionization and heating rates 
respectively due to the decaying dark matter particles. The additional terms include the 
contribution from both hydrogen (S = H) and helium (S = He) atoms. 

The hydrogen ionization rate Rm(z) due to the decaying dark matter particles is: 

R m (z) = nn(z) [1 - x R {z)\ f°° ^Pa H (E) dE, (4) 

where xh (z) is the ionization fraction at redshift z, %(z) is the total hydrogen number 
density including neutral and ionized hydrogen at redshift z, J (E) is the photon energy flux 
per unit energy per unit solid angle at energy E, and <th (E) is the photoionization cross 
section of hydrogen at energy E. The hydrogen heating rate Rhn(z) due to the decaying 
dark matter particles is: 

R m (z) = n H [l - *„(*)] ^ J W{E-E^) ^ {E) dK (5) 

J -Eth.H 

Similarly, the helium ionization rate Rm c (z) due to the decaying dark matter particles is: 



poo a t / rp\ 

Rme(z) = nn c (z) [1 - x Ue (z)] / — a Hc (E) dE, 

and the helium heating rate RhHe(z) due to the decaying dark matter particles is: 



(6) 



Rm e (z) = n Hc [l - x ne (z)} AnJ ^ E Eth ^ a Hc (E) dE. (7) 

The photon flux J(E) is 

m = _ [„ d „ o(1 + ,)3] , (8 ) 

where z em is the red shift when the dark matter particle decayed, E = E , 1 ^ rZ , E being 
the energy of the emitted photon, t(z) is the time elapsed since the big bang at redshift z, 
n dmo is the present number density of the decaying dark matter particles if they did not 
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decay so that ridmo(l + z ) 3 is the number density at redshift z, and r a b s is the optical depth 
caused by absorption of photons by ionization defined by 

1 + z' 

■)cdz' 

(9) 

(1 + z')#(0 



(1 + z')H(z') ^ ~ ^ ^ "h^W^ = 



+ 



2om | 1 + 2;' 

,\ H f A I 1 ~~ XHc CO) ^(^O^HeC-E = #0^ 



We consider photon energy E between the hydrogen threshold energy (=13.6 eV) 
and Eq at redshift z. The lower limit is £th,H since photons with E below -Eth.H do not 
have enough energy to ionize hydrogen. The upper limit is Eq since the redshift due to the 
expansion of the universe would cause photons to have energy E smaller than Eq. Hence 
the integral in (J3J) can be written as: 

AttJ(E) 



Eth.H 



E 



-cth (E) dE 



E 



c n dm0 (l + zf l— < 



^em)/^dm p T abs (^i^em) 



-a H (E) dE 



c n dm0 (l + z) 

7dm 



E a /(l+z) 



E tKn /(l+z) 



E H{z cm ) 

p —t(z eln =Eo / E ohs — l) /r dln p —r ahs (z,z eln =Eo / E obs — l) 



EnhnHiZp. 



Eq/E, 



obs 



-<r H (^obs(l + z))dE ohs , 
(10) 



where we have made a change of integration variable in the third line from E to E D ^ S , the 
present observed energy of a photon produced in the past by the decaying process. The 
integral for the helium contributi on can be rewritt en similarly. The photoionization cross 



section cr(E) is approximated by 



Osterbrockl ( 119741 ): 



<t(E) = a th lp 



E 



E 



th 



+ (1-/3) 



" E ' 




_Eth_ 





'ir 



where a th = 6.30 x 10~ 18 cm 2 , /3 = 1.34 and s = 2.99 for hydrogen, and cr th = 7.42 x 10~ 18 cm 2 , 
(3 = 1.66 and s = 2.05 for helium. To account for the energetic secondary electrons that 
could ionize and heat up hydrogen and helium atoms, we multiply t he cross sections in 
the ionization rates by an additional factor {1 + <f)[x(z)]E(z) / E^} (IMapelli &: Ferrara 
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20051 ). For the cross sections in the ionization rates, (f>(x) — C(l — x a ) b and C = 0.3908, 
a = 0.4092 and b = 1.7592 for hydrogen , and C = 0.0554, a = 0.4614 and b = 1.6660 



for helium (jShull &: van Steenberg 



19851 ). For the cross sections in the heating rates, 
(j){x) = C[l - (1 - x a f] and C = 0.9971, a = 0.2663 and b = 1.3163 for hydrogen and 
helium. In this model we take the Hubble parameter H(z) to be that given in the ACDM 
model, i?o\A^A + f2 m (l + z) 3 where Hq = H(z = 0) is the present Hubble parameter, Vt m is 
the present matter density, and f^A is the dark energy density. 

There are three free parameters md m , 7dm, / in this model, where md m and Td m are the 
mass and life time of the decaying dark matter particles, and / is the present mass fraction 
of the total dark matter accounted for by these hypothetical decaying particles. m dm and / 
are directly related to E and 7Zd m o m the above equations: 



E Q = m dm /2 

3gg 

idm 87rGm dm c 2 ' 



-f C~} ^o/^dn 



(12) 
(13) 



where Od m is the present dark matter density, G is the gravitational constant, and to is the 
age of the universe. If we assume Td m is much larger than to, the exponential factors in 
Eqs. (flOl) and ( 113]) can be ignored. The ionization and heating rates due to decaying dark 
matter particles then depend only on the ratio between / and r dm . 

By solving the two ordinary differential equations (0Q) and (J2J) we get the ionization 
fraction x and mat ter temperature T m as functions o f z. The subprogram RECFAST 



( jSeager et al. 



200(1 ) used in CAMB ( jLewis et al.l l2Q00h is modified to include the additional 



terms due to the decaying dark matter particles. We have changed CAMB so that the 
reionization history x(z) is calculated by our model, rather than the ad hoc tanh function 
built in. This introduces three new parameters: /, Td m and md m , which are varied together 
with other cosmological parameters. We have not changed the Boltzmann equations since 
the effects introduced by the decaying cold dark matter are very small and are of higher 
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order. For example, the photons from the dark matter decay affect other photons by first 
scattering with electrons which then scatter with other photons. Also, very few dark matter 
particles have decayed by the time of recombination. The fraction of electrons affected 
by the decay photons can be approximated by the ratio between the number densities of 
the decay photons and the CMB photons, which is at most of the order 1CT 7 within the 
interested ranges of / and r dm . Therefore only very few electrons and hen ce CMB photon s 
are a ffected by the decaying dark matter particles directly. COSMOMC (ILewis fc Bridle 



20021 ) and the 7-year WMAP data are used to constrain the parameters in this model. 



3. MCMC Results and Discussion 

By varying mam, 7dm , / and other standard parameters in COSMOMC and comparing 
the resulting CMBA spectrum with WMAP 7-year data, we obtain the constraints on these 
parameters. The lower limits of m dm , r dm and r dm / f are set to be 13.6 eV, 10 10 yr and 
10 10 yr respectively. The minimum value for Td m is chosen to be 10 Gyr because it has to be 
large enough so that not too much dark matter have already decayed by now. The result is 
shown in Figured! 

From the COSMOMC results the set of best fit parameters is shown in Table [U and 
the marginalized limits are shown in Table [2j 

We can see that m& m and r dm / / are highly correlated, and are less than about 1 MeV 



and 10 yr respecti vely. This mass 



models, for example 



Boyarsky et al. 



range 



overlaps with that of sterile neutrinos in some 



( 120061 ). For comparison, the ionization fraction using 
the set of best fit parameters in the decaying dark matter model and assumed in the original 
CAMB are plotted together in Figure |2l and the corresponding CMBA spectra are plotted 
in Figure |3j From Figure [2] we can see that the two different ionization fractions agree 
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23456789 
log 10 (m dm /eV) 



Fig. 1. — Left: Probability distribution of the parameters in the decaying dark matter model. 
The solid line is the marginalized likelihood and the dotted line is the mean likelihood. 
Right: Joint constraints of the three dark matter parameters in the decaying dark matter 
model. The contours refer to the marginalized likelihoods while the colours refer to the mean 
likelihood. The inner contour corresponds to 68% confidence limit and the outer contour 
corresponds to 95% confidence limit. Relative values of the likelihood are used, with the 
maximum likelihood set to be 1, and the color bar shows the relation between the colors 
and the relative likelihood. The 7- year WMAP results are used to constrain the parameters. 
Refer to Table [1] for explanation of the parameters. 9 is 100 times the ratio of the sound 
horizon to the angular diameter distance at about the time of recombination, and A$z is the 
floating amplitude for the SZ spectrum as in WMAP. 
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Table 1: Best fit parameters in the decaying dark matter model. A s is denned at the pivot 
scale of 0.05 Mpc" 1 . 



Parameter 


Symbol 


Value 


Hubble parameter (km s~ 1 Mpc^ 1 ) 


H 


69.6 


Physical baryon density 


Q b h 2 


0.02238 


Physical dark matter density 


Q c h 2 


0.1141 


Curvature fluctuation amplitude 


A s 


2.17 x 10~ 9 


Scalar spectral index 


n s 


0.963 


Decaying dark matter particle mass (keV) 




17.3 


Decaying dark matter particle life time over fraction (yr) 


Tdm/f 


2.03 x 10 16 



Table 2: List of the marginalized limits for different parameters in the decaying dark matter 
model, together with standard cosmological parameters. 



Symbol 


Prior 


Limits (68%) 


Limits (95%) 


H (km s^Mpc" 1 ) 


{40, 100} 


{67.5, 72.3} 


{65.3, 74.6} 


Q b h 2 


{0.005, 0.1} 


{0.02176, 0.02284} 


{0.02123, 0.02339} 


Q c h 2 


{0.01, 0.99} 


{0.1076, 0.1185} 


{0.1024, 0.1243} 


ln(10 10 A s ) 


{2.7, 4} 


{3.047, 3.115} 


{3.013, 3.149} 


n s 


{0.5, 1.5} 


{0.948, 0.975} 


{0.935, 0.987} 


log 10 (m dm /eV) 


{1.13, 9} 


{3.79, 4.84} 


{2.41, 6.05} 


log 10 (r dm ///Gyr) 


{1, 10} 


{6.10, 7.90} 


{3.36, 9.09} 
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quite well for most values of z, but the difference is significant for z ~ 50. Nevertheless 
the resulting CMBA spectra are still nearly the same. We can also see that the ionization 
fraction in the decaying dark matter model actually agrees quite well with the ad hoc 
imposed tanh function in default CAMB. Recently, joint CMBA-quasar absoprtion line 
constraints on the reionization history using a model independent principal component 
decomposition method suggests that reionization is 50 % complete betwe en 9.0 < z < 11.8, 



and 99% complete between 5.8 < z < 10.4 (95% CL) flMitra et al. 



2013 ). A similar study 



obtained a b est-fit reioniza t ion hi story very close to our result presented in Figure [2] 
(Figure 1 in 



Pandolfi et al 



(120111 )). Another recent work based on the patchy kinetic 
Sunyaev-Zel'dovich effect concludes that reionization ended at z > 5.8 or 7.2 (95% CL), 
depending on whe ther correlation with the cosmic infrared background is assumed or not 



(IZahn et al. 



20111 ). Our result shown in Figure [2] is consistent with these recent constraints. 



Our results are c onsistent with the constr aint from diffuse X-ra y background 



(IBoyarsky et al. 



20061 ). The empirical bound in IBoyarsky et al.l (120061 ) is log 10 (rd m /Gyr//) 



> 8.5, which gives further reduction of the allowed region in the contour plot in Figure [H 



4. Summary and Conclusion 

We have investigated the effects on CMBA by a component (mass fraction /) of dark 
matter particles with mass that decay with cosmological lifetime r^ m . The photons 
emitted are redshifted and may ionize hydrogen and helium at later times, affecting the 
reionization history of the universe. If Td m is much longer than the age of the universe, 
the optical depth depends only on the ratio of Td m and /. We obtained constraints on 
these parameters by using the WMAP 7-year data and modified RECFAST, CAMB and 
COSMOMC codes and assuming that the only reionization source is the decaying dark 



13 



matter. In the long lifetime limit, we find that 250 eV < m<i m < 1 MeV, 2.23 x 10 3 yr 
^ T dm/ f ^ 1-23 x 10 18 yr, and the best fit values of m dm and r dm / / are 17.3 keV and 
2.03 x 10 16 yr respectively. Sterile neutrinos with mass 17.4 keV are possible within our 



t 95% CL, wl 


lich r 


(Chan & Chu 


2011 



rich may account for the 8.7 



Prokhorov fc Silk 



teV emission observed by the 



20101 ) . The allowed range of Tdm/ f 



is reduced further if the constraint from diffuse X-ray backgro und is taken into account: 



3.16 x 10 17 yr < r dm // < 1.23 x 10 18 yr f lBovarskv et al. 



20061 ). 



We have shown that the reionization history of the universe is sensitive to decaying 
dark matter parameters, and future experiments may lead to tighter constraints on dark 
matter models. 



This work is partially supported by grants from the Research Grant Council of the 
Hong Kong Special Administrative Region, China (Project Nos. 400805 and 400910). We 
thank the ITSC of the Chinese University of Hong Kong for providing its clusters for 
computations. 
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Fig. 2. — Left: The upper figure plots the ionization fraction x(z) versus z, and the lower 
figure plots the matter temperature T m {z) versus z, both using the decaying dark matter 
model with the best fit parameters in Table [U and the original CAMB with the best fit 
WMAP parameters. Right: Same as left frame, but for the region around z ~ 10. 
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Fig. 3. — CMBA temperature power spectra, calculated using the decaying dark matter 
model with the best fit parameters in Table [U and the original CAMB with the best fit 
WMAP parameters. The values of the best fit — ln(likelihood) in the decaying dark matter 
model and standard ACDM model are 5527.44 and 5532.39 respectively. 
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